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The temperature dependence of the uniform susceptibility 
and the ground state energy of antiferromagnetic Heisenberg 
ladders with up to 6 legs has been calculated, using the Monte 
Carlo loop algorithm. The susceptibilities of even-leg-ladders 
show spin gaps while these of odd-leg-ladders remain finite in 
the zero temperature limit. For small ratios of intra- to inter- 
leg couplings, odd-leg-ladders can be mapped at low temper- 
atures to single chains. For equal couplings, the logarithmic 
corrections at low temperatures increase markedly with the 
number of legs. 

Recently, antiferromagnetic Heisenberg spin-1/2 lad- 
der systems have attracted much interest, following the 
discovery of a finite spin gap in the 2-leg-ladder . Later 
investigations showed that the crossover from the sin- 
gle Heisenberg chain to the two-dimensional (2D) an- 
tiferromagnetic square lattice, obtained by assembling 
chains to form "ladders" of increasing width, is far from 
smooth 1^. Heisenberg ladders with an even number of 
legs (chains), ni, show a completely different behavior 
than odd-leg- ladders. While even- leg-ladders have a spin 
gap and short range correlations, odd- leg-ladders have no 
gap and power-law correlations. Based on density matrix 
renormalization group (DMRG) studies. White et al. |] 
gave an explanation of this fundamental difference in the 
framework of the Resonant Valence Bond (RVB) picture 
These theoretical predictions have been verified ex- 

and the 
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perimentally, in materials such as {yO)2P20j \ 
homologous series of cuprates Srn-iCun+i02n 
contain weakly coupled arrays of ladders. 

Previous numerical studies, using exact diagonaliza- 
tion the Quantum Monte Carlo world line algo- 

rithm |g[ or the Quantum Transfer Matrix method 
were restricted to small systems or could not be applied 
at low temperatures. Using the new Quantum Monte 
Carlo (QMC) loop algorithm we overcome these lim- 
itations and are able to investigate very large ladders to 
much lower temperatures. 

In this letter we consider ladders with ni legs (n; = 
1,2,..., 6) of length L. The Hamiltonian of such 
spin— 1/2 systems 
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is defined on L x ni lattices. The sum marked by <-> (|) 
runs over nearest neighbors along the legs (rungs). We 
assumed periodic boundary conditions in the longitudi- 
nal direction of the ladder. The exchange constants J 



and J± are positive, corresponding to antiferromagnetic 
coupling. 

Using the Quantum Monte Carlo loop algorithm with 
improved estimators, we have calculated the tempera- 
ture dependence of the uniform susceptibility x a-nd the 
internal energy E. The QMC loop algorithm was first 
developed by Evertz et al. for the 6-vertex-modcl, 
but can also be applied to quantum spin systems pd| , p^ . 
The QMC loop algorithm is an improved world line algo- 
rithm. The updates in the loop algorithm are global and 
no longer local as in the conventional Metropolis world 
line algorithm. This has the great advantage that the 
autocorrelation times are reduced by several orders of 
magnitude. It was thus possible to simulate very long 
ladders to very low temperatures. We considered sys- 
tems up to 100 x 6 sites and reached temperatures down 
to T = J/50 without major problems. All results are 
extrapolated to a Trotter time interval At 0. The 
application of improved estimators (see e.g. Ref. |^), 
further reduces the variance of the measured observables 
dramatically. 

In Table | the ground state energies of the different 
ladders in the isotropic case (J — J±) are presented. 
Considering ladders of lengths L, we first extrapolate to 
T ^ 0. For the odd-leg-ladders we use the form 
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where E^ (T) is the internal energy of the ladder of length 
L. This form is motivated by the infinite single chain, 
which in a low temperature field theory can be described 
by a massless boson. Our numerical data agree well with 
Eq. (^) for the finite single chains, as well as for the fi- 
nite 3- and 5-leg-ladders. The even-leg-ladders, on the 
other hand, have spin gaps. The internal energy for T 
well below the spin gap A is determined by the thermal 
occupation of the lowest lying S = 1 magnon band with a 
quadratic dispersion near the zone boundary minimum. 
For the extrapolation to T — > we therfore use the form 
ElIT) = El{Q) + b (T3/2 + 2AT1/2) exp(-A/T). In a 
second step the ground state energies El{Q) for the finite 
systems are extrapolated to the bulk limit i?L=oo(0), fit- 
ting El{Q) to a polynomial in 1/L. The finite size correc- 
tions of the ground state energy, however, are negligibly 
small for ladders with L > 100. 

In Fig. ^ we show the susceptibility per rung of ni 
spins, xi'P)- for fti6 ladders with 1 < < 6 in the 
isotropic case J = J±. We always consider ladders long 
enough such that finite size corrections are negligible. For 
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T > J the results agree well with a third order high tem- 
perature expansion. At low temperatures, we observe 
the predicted behavior. Ladders with even ni show an 
exponential drop of the susceptibility indicating a gap 
in the excitation spectrum. For larger nj, the drop sets 
in at smaller T and is steeper. The gap A„j decreases 
substantially with increasing ni . For odd n; on the other 
hand, x(r) remains finite also at T ^ J, as in the single 
chain. 

Ladders with even ni were already investigated in de- 
tail PJs|Jl5[|. For temperatures T <^ A the susceptibility 
for the 2-leg-ladder is determined by the thermal occupa- 
tion of S" = 1 magnon band with a quadratic dispersion 
near the zone boundary minimum M : 



X{T) cx T 
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Provided a quadratic dispersion for the lowest lying 
magnon branch in the excitation spectrum near its mini- 
mum is assumed, Eq. (||) also holds for the 4- and 6-leg- 
ladder. We estimate A„, by fitting the numerical QMC 
data for low temperatures and find in the isotropic case 
A2 = 0.51(1)J for ni = 2, A4 = 0.17(1)J for n/ = 4 
and Ae = 0.05(1)J for ni ~ 6. The value A2 obtained 
for the 2-leg-ladder is in perfect agreement with former 
results 1^,^,^. On the other hand, the spin gap obtained 
by White et al. §, using DMRG methods for the 4-leg- 
ladder A4 = 0.190 J is slightly larger than our value. 

The decrease of the spin gap with increasing ni can 
be explained by delocalization of RVB singlets not only 
along but more and more also across the ladder. The 
decrease of the spin gap, however, is much faster than 
A„, cx 1/ni, suggested in The spin gap for the 6- 
leg-ladder Ag is already a factor 10 smaller than A2, 
suggesting rather an exponential decrease of A„, with 
increasing n;. 

The susceptibility per rung of the odd-leg-ladders re- 
mains finite in the low temperature limit and tends to a 
zero temperature value approximately independent of ni 
see Fig. 11(b)]. For this we conclude that the odd-leg- 
ladders belong to the same universality class as the single 
chain. 

The single chain can be described in a low temperature 
field theory by the k = 1 Wess-Zumino-Witten nonlinear 
cr-model or equivalently by a free, massless boson, with a 
(spin) velocity v = 'kJ/2. Based on this model x{T = 0) 
and the leading T dependences of x have been calculated 
p6| , p7t with the leading marginally irrelevant operator. 
Two-loop renormalization of the marginal coupling leads 
to 



ln(To/T) 

+ 0((ln T)-3) 



ln(ln(To/T)-H/2) 
21n^(To/T) 



(4) 



where Tq is the cut-off-temperature. The susceptibil- 
ity approaches its asymptotic zero temperature value 



X(0) = (27ri;)-i = (Jtt^)"! with infinite slope. The field 
theoretical results can be compared to the exact Bethe 
ansatz data and one finds that Eq. (^ holds to 
within 1% for T < O.IJ. These results are shown in Fig. 
H together with our QMC data for low temperatures. 

We turn now to the 3- and 5-leg-ladders. In the limit 
J/ J±_ — each eigenfunction is a direct product of one- 
rung states whose lowest lying multiplet is a spin doublet 
(s = 1/2). The ground state of the whole system is 
therefore 2^-fold degenerate. A finite value of J lifts this 
degeneracy. In this 2^-dimensional subspace Ai we can 
define an effective Hamiltonian Hes which includes all 
intra-leg interactions. To first order in J/J± we obtain 
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Sj. totSj+l, tot, 
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where Sj^ tot is the total spin of the j rung and Joff = J 
for the 3-leg-ladder, respectively Joff = 1.017 J for the 
5-leg-ladder. has just the form of the Hamiltonian 

of the single chain with an effective coupling Joff and we 
can map the low lying energy states of the 3- and 5-leg- 
ladder to those of the single chain. In the following we 
will concentrate on the 3-leg-ladder. 

The susceptibility of the single chain xi(T/J) scales 
with 1/J. It follows, that for a 3-leg-ladder with small 
J / J±_ and at low temperature, where only the above men- 
tioned low lying states in M. are relevant, the suscepti- 
bility per rung xs scales with 1/Jcff and has the same 
functional dependence on T/Jcff as xi{T/J), according 
to Eq. (|). This can be seen in Fig. ^ where we show 
JX3 for different ratios J / J±_ as a function of T/ J. For 
small J / J± the susceptibility per rung xs multiplied by J 
is very close to Jxi until a crossover temperature, which 
depends on J/J±. Above this temperature, the suscep- 
tibility of the 3-leg-ladder is larger, due to the presence 
of additional states in the 3-leg-ladder which are not in- 
cluded in the 2^-dimensional subspace Ai. These addi- 
tional states have a finite gap A. The susceptibility of 
the 3-leg-ladder then reads 



X3 = Xl(^cff) + X, 



(6) 



where x is the contribution of the additional states. From 
our QMC data we find JcS ~ J for all small J/J± and 
A4(J/J±) < A < A2(J/Ji). 

In the isotropic case J ~ J±, the gap A in units of J 
is indeed smaller than for small J/J±_ but remains finite 
(A4 < A < A2). For r <C A we can neglect the contri- 
bution X of the additional states. Comparing X3 to xi 
for T ^ A, we see that their slopes are completely differ- 
ent (see Fig. ^). Therefore we conclude that in this case 
the simple model, discussed above [Eq. (^], no longer 
applies. 
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As J/Jl — > 1 also next-nearest neighbor and longer 
range interactions between rung-spins become important 
in the effective Hamiltonian fleff- Since these additional 
interactions respect the SU (2) and translational symme- 
try, Eq. (||) still applies for some values of v and Tq, 
according to p6[ . Therefore, for very low temperatures 
the susceptibility of the 3-leg-ladder can be described by 
Eq. (Q) also in the isotropic case (see Fig. ^). 

The spin velocity in the isotropic 3-leg-ladder seems 
to be close to that of the single chain vi . This can be seen 
by two facts. First, both of the susceptibilities per rung 
seem to extrapolate to the same zero-temperature value. 
Secondly, White et al. Q determined the spin gap of the 
finite single chain and the finite 3-leg-ladder in function 
of 1/L. The slopes of these curves as 1/L — > 0, ttvi and 
7rw3, agree, at least within 5%. Assuming vi — W3, we 
get a rough estimate of the cut-off-temperature Tq in the 
isotropic 3-leg-ladder. The value is much smaller than 
in the single chain (see Fig. ||). We conclude therefore, 
that the effective interactions between the spinous in the 
3-leg-ladder are much stronger than those in the single 
chain. 

We conclude that the odd-leg-ladders belong to the 
same universality class as the single chain and can be 
described in the zero temperature limit by a k=l Wess- 
Zumino-Witten non linear cr-model with a spin velocity 
Vni ■ These velocities seem to have more or less the same 
value for all n/. With increasing n;, however, we move 
further away from the conformal point and the logarith- 
mic corrections, due to the leading marginally irrelevant 
operator, increase markedly. 

Finally, we want to point out, that the zero tempera- 
ture value of the susceptibility per site X"/*''^ (0) for the 
odd-leg-ladders decreases with increasing ni. Since the 
odd-leg-ladders seem to have more or less the same zero 
temperature value of the susceptibility per rung, it fol- 
lows that Xni'*°^(0) For — > 00 the zero tem- 
perature value Xn (0) therefore goes to zero for odd ni 
as well as for even nj. The susceptibility per site of a 
2D square lattice, however, is finite for T = 0. This is 
therefore a further example that the crossover from the 
single chain to the 2Z?— lattice is not a smooth one. 
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number of legs obtained ground reference value 

state energy 

1 -0.4432(1) -0.44315... 

2 -0.5780(2) -0.578 §] 

3 -0.6006(3) 

4 -0.6187(3) 

5 -0.6278(4) 

6 -0.635(1) 

2D lattice -0.6693(1) 0| 



TABLE I. Ground state energies per site for the differ- 
ent ladders in the isotropic case. For the single chain we 
have perfect agreement with the analytical result | — In 2 
from the Bethe Ansatz ^]. Furthermore, the result for the 
2-leg-ladder coincides with the ground state energy calculated 
by Barnes et al. using Lanczos techniques. With increas- 
ing width the results approach to the ground state energy per 
site of the infinite 2D square lattice, which was calculated 
by various methods. For an overview see jl^. The reference 
value given here was recently obtained by U.J. Wiese and H.P. 
Ying ||l2| using the QMC loop algorithm. 



FIG. 1. Susceptibility as a function of the temperature for 
the single chain and the Heisenberg ladders with up to 6 legs. 
At high temperatures the result agree well with a third order 
high temperature expansion. The low temperature region is 
shown in (b) in a larger scale. To distinguish the curves some 
data points are marked by symbols. The error bars are smaller 
than the symbols. 



FIG. 2. Renormalization group improved field theory (solid 
lines) [Eq. (^)] for different cut-off-temperatures To versus 
Bethe ansatz data and QMC results for x{T) at low tem- 
perature. The error bars are smaller than the symbols. 



FIG. 3. Susceptibility per rung of the 3-leg-ladder xs 
(dashed lines) for different ratios J/ J± and of the single chain 
Xi in function of the temperature. The inset shows the low 
temperature region in a larger scale. To distinguish the curves 
some data points are marked by symbols. The error bars are 
smaller than the symbols. 
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